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Abstract 

We have defined pinning fields as those random fields that keep some of the magnetic 
moments unreversed in the region of negative external applied field during the 
demagnetizing process. An analysis of the statistical properties of such pinning 
fields is presented within the context of the Gaussian Random Field Ising Model 
(RFIM). We show that the average of the pinning fields exhibits a drastic increase 
close to the coercive field and that such an increase is discontinuous for low degrees 
of disorder. This behaviour can be described with standard finite size scaling (FSS) 
assumptions. Furthermore, we also show that the pinning fields corresponding to 
states close to coercivity exhibit strong statistical correlations. 
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1 Introduction 



The magnetization reversal process in a ferromagnet is a complex dynamical 
process which is still not totally understood [1] . It is, nevertheless, very impor- 
tant for both fundamental and technological reasons. Different factors play a 
role in determining the metastable path that the system follows to reverse the 
magnetization from full positive saturation to full negative saturation when 
the external field is decreased. Among others, thermal fluctuations, long range 
dipolar forces, anisotropy and local forces due to disorder, compete together 
in order to decide the sequence of magnetic domains that transform. 

A first simplification of the problem consist in neglecting the role of fluctu- 
ations and relaxation effects. This corresponds to the limiting case of "rate- 
independent" hysteresis. Magnetization reversal steps occur as almost instan- 
taneous avalanches joining metastable states. Such avalanches are triggered 
when the external forces induced by the applied field are strong enough to 
overcome the internal energy barriers which are caused by exchange forces, 
long-range dipolar forces and forces created by disorder. 

A prototypical model for the study of the influence of disorder in such an 
"athermal case" is the Gaussian RFIM. This model has been studied using 
two different approaches. On the one hand, a number of studies[2,3] have fo- 
cused on the analysis of a single magnetic interface. The numerical algorithms 
assume that only spins close to the interface can flip. On the other hand, stud- 
ies including both nucleation events and interface motion have been performed 
by using synchronous relaxation dynamics [4,5]. In both cases hysteresis ap- 
pears as a consequence of the local fields that keep the magnetic moments 
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unreversed (pinned) even at negative values of the applied external field. We 
will focus on the statistical analysis of such pinning fields in the case of RFIM 
with synchronous relaxation dynamics. 

Our goal is to point out an essential difference between the two approaches. 
The quenched pinning fields originating from disorder exhibit very different 
statistical distributions in both cases. The distribution of pinning fields in an 
intermediate state is very different from the quenched disorder distribution 
corresponding to the initial saturated state (Gaussian) in lattice models that 
include nucleation and interface movement. This is due to the fact that in 
the initial stages of the demagnetization process the regions with low energy 
barriers have already been reverted. Close to coercivity, the remaining barriers 
are much higher due to the previous selection process. Such an effect does not 
occur for the models with an advancing single interface, for which the pinning 
fields in the unreversed regions always exhibit a statistical distribution which 
corresponds to the original quenched disorder distribution. 



2 Model 



The 3d-Gaussian RFIM at T = is defined on a cubic lattice of size N = L 3 . 
On each lattice site a variable Si — ±1 accounts for the magnetic degrees of 
freedom. The Hamiltonian is: 

H = -'£S i S j - y £S i h i -B y £S i (1) 

ij i i 

where the first term stands for the exchange interaction between nearest- 
neighbour spins Si, the second term for the interaction with the quenched 
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local random fields hi and the last term for the interaction with the driving 
field B. The hi are Gaussian distributed with zero mean and variance a 2 . 
Although this model enable long-range dipolar forces to be included, from 
a computational point of view the numerical solution becomes much harder. 
Since we are interested in the analysis of the pinning fields generated by dis- 
order, we have neglected long-range terms. 

The numerical simulations are performed using local relaxation dynamics [4]. 
The initial saturated state with all the spins Si — 1 corresponds to the equi- 
librium state with B = +00. The field B is then decreased until a spin 5, 
becomes locally unstable. At this point the external field is kept constant and 
the spin Si is flipped. This may cause an avalanche since some of its neigbour- 
ing spins may become unstable. All unstable spins are flipped synchronously, 
until the avalanche ends. The external field B is then decreased again until a 
new avalanche starts. 

The hysteresis loop is obtained by measuring the magnetization m = X^Li Si/N 
as a function of B. Fig. 1 shows an example and the corresponding configu- 
rations snapshots. As can be seen, nucleation and interface movement coexist 
during evolution. Due to the finite size of the simulated system, the loops 
consist of a sequence of discontinuous jumps or avalanches for each realiza- 
tion of the random fields hi. As has been studied in previous works [6,7], the 
characteristics of the loops depend on the amount of disorder a in the ther- 
modynamic limit (L — > 00): they are continuous for a > a c , whereas they 
display a discontinuity (corresponding to a spanning avalanche) at the coer- 
cive field for a < a c . This behaviour is associated with the existence of a 
"metastable" critical point on the (B, a) phase diagram located at a c = 2.21 
and B c = —1.425. The behaviour close to this critical point can be described 
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by a set of critical exponents. For instance, the correlation length diverges 
with an exponent v ~ 1.2, the order parameter (the magnetization jump Am) 
goes to zero with an exponent (5 ~ 0.024 when a — > a c from below and, at 
a = a c , the magnetization behaves as m ~ \B — B^ 1 ^ 5 with 5 ~ 50. Such crit- 
ical exponents have been obtained by detailed FSS analysis[6] which has also 
revealed that the most convenient scaling variable that measures the distance 



to a c is 



U U C A 

u = + A 



with A = —0.2. Furthermore, since we will be interested in the measurement 
of properties as a function of the external field B, we will need a second scaling 
variable to measure the distance to B c . The first simpler choice is: 

B-B r 



-Dcoe 



where B coe (cr,L) is the coercive field that tends to B c when a — > a c and 
L — > oo. 



3 Results 

We define pinning fields hf as those quenched random fields hi for which 
Si = +1 during the reversal process for a certain intermediate configuration. 
Although the set of pinning fields is simply a subset of the original quenched 
random fields, their statistical properties depend on the exact path followed 
until a certain configuration is reached. We have computed the average value 
of the pinning fields (hf) and the histograms corresponding to their statistical 
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distribution f(hf) as a function of B and a by simulating many realizations 
of disorder. Moreover we have measured pair correlations as C(hf,hf) = 
L*((hth+)-(ht)*). 

Fig. 2 shows the evolution of the distribution of pinning fields along the de- 
creasing branch of a hystereis loop corresponding to a = 2.14 < a c . This 
distribution can be understood as the distribution of barriers created by the 
quenched disorder that keeps the spins in the metastable state. Although ini- 
tially the distribution of pinning fields is similar to the original Gaussian dis- 
tribution, as the magnetization decreases, f(hf) starts to develop a non-trivial 
structure. In general the distribution tends to shift to the right, towards the 
region of large pinning fields, but it also develops a number of peaks associated 
with the seven possible local magnetization environments. 

Figs. 3 and 4 display the evolution of (hf) and C(hf,hf) as a function of 
the external field B for two different values of a corresponding to two cases: 
above and below a c . It is interesting to point out that the correlation is not 
zero in the two cases and displays a peak close to the coercive field. Note 
that this means that the distributions f(hf) are only projections of complex 
multivariate distributions. Additionally, (hf) displays a discontinuity for a < 
a c . The behaviour of (hf) is, therefore, similar to the behaviour of an order 
parameter around a critical point. 

Quantitative analysis of such data requires a convenient FSS analysis to be 
carried out. The two analyzed properties (hf) and C(hf, hf) are functions of 
the external field B, the amount of disorder a and the system size L. The FSS 
ansatz allows to the singular (critical) contributions of these two properties to 
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be expressed as a function of the invariants x = uL x l u and y = vL^l u : 



The exponents 9 and 3 — p characterize how (hf) decreases and how the 
correlation C diverges. Figures 5 and 6 show the scaling functions h and C. 
The good quality of the collapses of data corresponding to different system 
sizes, confirms the scaling assumptions. The exponents that allow the best 
collapses are 9 ~ 0.41 and 3 — p ~ 1.80. Moreover, the behaviour of the 
scaling functions allows prediction of how the average pinning field and the 
correlations behave in the thermodynamic limit. Since h behaves as -\f v l^ for 
v < (as indicated by the discontinuous line) , (hf) is finite for B < B c . 
The scaling behaviour for v > is not so good since there are non-scaling 
contributions due to the existence of non-spanning avalanches for all values of 
a [6] . As regards correlations, the scaling functions display a peak profile which 
indicates that, besides non-scaling contributions, the correlation C{hf ,h1~) 
will diverge at B = B c and a = cr c in the thermodynamic limit. 

4 Conclusions 

Pinning fields are responsible for the energy barriers that keep the spins in the 
metastable state within the context of the 3d RFIM at T = 0. The statistical 
distribution of pinning fields during the reversal magnetization process has 
been studied. Initially, when the system is saturated, the pinning fields are 
trivially distributed according to the nominal Gaussian distribution of random 
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fields. As the demagnetizing process advances and domains of negative spins 
are created, the pinning fields display a complex distribution. Their mean 
value increases monotonously for decreasing B. For a > a c this increase is 
continuous but for er < a c the average pinning field displays a discontinuity at 
coercivity. Moreover we have shown that close to coercivity the pinning fields 
exhibit strong statistical correlations. We finally remark that such a complex 
behaviour of the distribution of pinning fields is not taken into account in the 
studies that focus on the analysis of an advancing single magnetic interface. 

We acknowledge fruitful discussions with F.J.Perez-Reche. This work has re- 
ceived financial support from CICyT (Project No MAT2001-3251) and CIPJT 
(Project 2001SGR00066). X.I. aknowledges financial support from DGI-MCyT 
(Spain). 
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Fig. 1. Example of a hysteresis loop corresponding to a = 2.5 and L = 60. The 
insets show examples of typical spin configurations. 
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Fig. 2. Examples of distributions f(hf) corresponding to L = 60, a = 2.14 and 
different values of the external field as indicated. Data corresponds to averages over 
1000 realizations of disorder. 
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Fig. 3. Example of the evolution of (hf) and C(hf, h^) as a function of the external 
field for a system with L = 60 and a = 2.28. Data correspond to averages over 1000 
realizations of disorder. 
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Fig. 4. Example of the evolution of (hf) and C(hf, h^) as a function of the external 
field for a system with L = 60 and a = 2.14. Data correspond to averages over 1000 
realizations of disorder. 
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Fig. 5. FSS analysis of the average pinning field (hf). The figures show, on log-log 
scales, the scaling function h in front of y = vL^ 8 l v for 4 different values of 
x = ul}l v as indicated on each graph. Data correspond to the overlap of sizes 
L = 15, 30, 60, 120 (as indicated by the legend) and averages over many realizations 
of disorder ranging from 10 6 to 10 3 . 
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. 6. FSS analysis of the correlation between pinning fields C(hf, hj~). Details 
same as in Fig. 5. 
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